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Princeton  University,  Princeton,  NJ  08544 
609-258-5271  (Voice);  609-258-6233  (Fax);  cklaw@princeton.edu 

Statement  of  Problem 

The  present  program  aimed  to  study  ignition  in  nonhomogeneous  systems  that  are  of 
relevance  to  practical  combustion  devices.  The  investigation  was  conducted  in  both 
laminar  and  turbulent  environments,  and  involved  experimental,  computational,  and 
analytical  components. 

For  turbulent  ignition,  four  specific  projects  were  xmdertaking,  namely:  (1) 
Characterization  of  non-reacting  turbulent  flow  fields  in  counterflow,  in  which  the 
velocity  moment  along  the  counterflow  axis,  mean  and  fluctuating  velocity  profiles,  and 
turbulent  power  spectrum  were  measured  for  a  wide  range  of  pressures,  bulk  strain  rates, 
perforated  plate  configurations,  and  air  temperatures;  (2)  turbulent  non-premixed  ignition 
experiments,  for  which  ignition  temperatures  were  determined  for  hydrogen  for  a  range 
of  turbulent  flows,  fuel  concentrations,  and  pressures.  It  was  observed  that  optimal 
conditions  for  ignition  exist  in  weakly  turbulent  flows  where  the  ignition  temperature  is 
lower  than  either  laminar  or  strongly  turbulent  flows  at  the  same  pressure;  (3)  modeling 
study  of  turbulent  nonpremixed  ignition,  which  was  conducted  by  using  Monte  Carlo 
techniques  to  solve  a  joint  scalar  PDF  equation  to  assist  the  rmderstanding  of  our 
turbulent  ignition  experimental  results;  (4)  computational  simulation  of  the  ignition  of  a 
nonpremixed  hydrogen/air  mixing  layer  with  an  embedded  vortex  with  detail  chemistry, 
showing  that  the  scalar  dissipation  rate  is  a  crucial  controlling  factor  for  ignition,  causing 
a  nonmonotonic  dependence  of  ignition  delay  on  various  vortex  parameters. 

In  parallel,  three  projects  in  laminar  premixed  and  nonpremixed  coimterflows 
were  explored,  especially  focusing  on  the  chemical  kinetics  aspects  and  the  coupling  with 
transport,  namely:  (1)  Premixed  hydrogen/air  and  propane/air  ignition,  showing  that 
diffusion  can  modify  the  local  ignition  kernel  conditions  to  be  different  with  the 
boundary  conditions,  and  as  a  consequence  the  hydrogen  premixed  ignition  exhibits  five 
limits  in  the  pressure-temperature  diagram;  (2)  experimental  determination  and 
computational  calculation  of  the  counterflow  ignition  temperatures  and  laminar  flame 
speeds  of  several  C2-C3  hydrocarbons,  showing  deficiencies  in  literature  mechanisms;  (3) 
experimental  determination  and  computational  calculation  of  the  ignition  temperatures  of 
nonpremixed  counterflow  dimethyl  ether  (DME)  versus  air,  leading  to  the  validation  of  a 
recent  mechanism  that  showed  very  close  agreement  with  the  experimental  data. 
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Summary  of  Results 


1.  Effects  of  Turbulence  on  Ignition 

1.1  Characterization  of  non-reacting  turbulent  flow  fields 

High  fidelity  experimental  data  for  turbulent  non-reacting  counterflows  are  necessary  to 
interpret  turbulent  non-premixed  ignition  experiments  and  to  develop  and  evaluate 
numerical  models  for  turbulent  ignition.  It  is  particularly  important  to  have 
measurements  of  counterflows  where  one  jet  has  a  significantly  different  density  than  the 
other,  as  such  data  is  not  available  in  the  literature.  Counterflowing  jets  of  heated  air  and 
cold  nitrogen  were  investigated  using  a  two-component  laser  Doppler  anemometry 
system  to  measure  velocity  moments,  spectra,  and  autocorrelations  throughout  the  flow 
field.  Operating  conditions  spanned  a  range  of  pressures,  bulk  strain  rates,  perforated 
plate  configurations,  and  air  temperatures.  In  addition  to  measuring  velocity  moments 
along  the  counterflow  axis,  mean  and  fluctuating  velocity  profiles  were  made  in  the  radial 
direction  at  each  nozzle  exit.  These  profiles  are  important  for  establishing  boundary 
conditions  in  turbulent  flow  models.  The  radial  velocity  gradients  are  needed  in  quasi- 
one  dimensional  turbulence  models  and  the  profiles  themselves  are  needed  for  higher 
dimensional  Reynolds  stress  models  and  Large  Eddy  Simulations. 

Autocorrelation  measurements  at  the  burner  exits  were  used  with  Taylor's 
hypothesis  to  estimate  the  flow's  integral  length  scale.  The  length  scales  were  about  70% 
of  the  perforated  plate  hole  diameters,  in  agreement  with  theory  and  other  experiments. 
The  length  scale  is  important  for  modeling  efforts  since  it  can  be  used  to  estimate  the 
turbulent  dissipation  rate. 

Turbulent  power  spectra  were  determined  at  the  burner  exits  to  fiirther 
characterize  the  turbulence  for  a  given  perforated  plate  configuration.  For  example,  due 
to  the  difficulties  associated  with  generating  turbulence  in  a  strongly  heated  flow,  some 
experiments  have  been  conducted  with  only  one  perforated  plate,  located  in  the  cold 
flow.  The  use  of  this  configuration  was  justified  by  the  observation  that  the  turbulent 
fluctuations  resemble  those  in  a  turbulent  counterflow  that  has  perforated  plates  in  both 
nozzles.  However,  a  comparison  of  the  power  spectra  measured  at  the  hot  side  of  the 
burner  shows  that  there  is  a  significant  difference  between  the  two  cases.  The  spectrum 
of  the  case  without  a  second  perforated  plate  is  suggestive  of  very  young,  undeveloped 
turbulence,  indicating  that  it  is  important  to  have  a  perforated  plate  in  both  sides  of  the 
burner  when  studying  turbulence. 

1.2  Turbulent  non-premixed  ignition  experiments 

Our  previous  studies  on  turbulent  nonpremixed  ignition  were  performed  in  the 
counterflow  configuration,  with  turbulence  generated  by  a  perforated  plate  in  the  fuel 
side  of  the  flow  as  was  done  by  other  studies  involving  counterflow  turbulent 
combustion.  For  the  hydrogen/air  system,  it  was  found  that  the  turbulence  intensity  of 
the  flow  did  not  have  an  appreciable  effect  on  the  ignition  temperature.  These  results 
were  understood  to  reflect  the  facts  that,  in  laminar  flows,  hydrogen  ignition  is  relatively 
insensitive  to  strain  rate  variations  under  second-limit  conditions  such  as  those  present  at 
atmospheric  pressure,  and  that  the  role  of  turbulence  in  the  flows  studied  was  felt 
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primarily  in  the  form  of  a  turbulent  strain  rate.  However,  since  in  these  experiments 
turbulence  was  generated  only  on  the  fuel  side,  while  the  ignition  kernel  was  located  on 
the  oxidizer  side,  it  is  not  certain  how  representative  the  results  are  of  real  turbulent 
conditions.  Thus  a  major  objective  of  the  present  investigation  was  to  study  ignition  with 
turbulence  generated  on  both  sides  of  the  flow  by  using  two  perforated  plates.  Results 
that  are  qualitatively  different  from  the  previous  work  were  observed  and  are  reported  in 
the  following. 

Ignition  temperatures  were  determined  for  hydrogen  for  a  range  of  turbulent 
flows,  fuel  concentrations,  and  pressures.  In  the  pseudo-turbulent  configuration  of  one 
perforated  plate,  the  flows  are  characterized  by  a  narrow  range  of  large  scale,  low 
frequency  fluctuations  that  produce  a  quasi-steady  response.  With  perforated  plates  on 
both  sides  of  the  flow,  the  lowest  turbulence  intensity  investigated  significantly  reduced 
ignition  temperatures  at  atmospheric  pressure  relative  to  laminar  values  as  shown  in  Fig. 
1.  Subsequent  increases  in  turbulence  intensity  caused  the  ignition  temperatures  to  rise 
and  the  dependence  on  the  bulk  strain  rate  to  weaken.  These  results  indicate  that  optimal 
conditions  for  ignition  exist  in  weakly  turbulent  flows  and  these  temperatures  are  notably 
lower  than  can  be  obtained  at  the  same  pressure  in  laminar  flows.  Compared  to  pseudo- 
turbulent  flows,  fully  turbulent  flows  have  a  wider  range  of  eddy  sizes  that  extend  to 
much  smaller  scales,  thereby  affecting  the  structure  of  the  ignition  kernel  and  the  ignition 
behavior. 

The  effects  of  varying  the  fuel  concentration  in  turbulent  flows  were  found  to  be 
qualitatively  similar  to  those  in  laminar  flows;  ignition  is  relatively  insensitive  to  fuel 
concentration  at  high  fuel  concentrations  and  becomes  increasingly  more  sensitive  at 
lower  concentrations.  The  relative  order  of  the  ignition  temperatures  between  different 
flows  was  the  same  for  all  fuel  concentrations. 

The  effect  of  pressure  on  ignition  is  of  particular  interest  due  to  the  prominently 
nonlinear  nature  of  the  hydrogen  oxidation  kinetics  as  exemplified  by  the  p-T  explosion 
limits.  Further  experiments  were  conducted  at  three  turbulence  intensities  and  pseudo- 
turbulent  flows  at  pressures  ranging  from  0.6atm  to  8atm  and  a  fuel  concentration  fixed  at 
20%.  Results  are  shown  in  Fig.  2  for  two  different  pressure-weighted  strain  rate  ranges, 
with  laminar  calculations  and  crossover  temperatures  included  for  comparison.  The 
increase  in  ignition  temperature  associated  with  increasing  hrrbulence  intensity,  noted 
above  for  one  atmosphere,  is  seen  to  hold  for  all  pressures  studied.  Furthermore,  the  large 
decrease  in  ignition  temperature  from  laminar  ignition  to  the  lowest  turbulence  intensity 
is  found  to  hold  for  pressures  throughout  the  second  ignition  limit.  However,  at  higher 
pressures,  the  laminar  ignition  temperature  increases  more  slowly  as  the  transition  occurs 
between  the  second  and  the  third  ignition  limits  so  that  at  a  fixed  pressure  the  turbulent 
ignition  temperature  becomes  higher  than  the  laminar  calculation.  At  the  lower  bulk 
strain  rates  of  Fig.  2a,  the  turbulent  ignition  response  seems  to  undergo  a  transition  to  the 
third  limit  where  the  ignition  temperature  decreases  with  increasing  pressure,  although 
this  transition  occurs  at  6atm.  in  the  turbulent  flows  rather  than  at  5atm.  in  a  laminar  flow 
with  a  similar  strain  rate.  The  data  at  the  higher  bulk  strain  rates  in  Fig.  2b  do  show  the 
tendency  to  transition  to  the  third  limit,  although  without  a  clear  turning  in  the  trend.  The 
observed  behavior  of  turbulent  flows  represents  an  extension  of  the  second  limit  behavior 
to  higher  pressures  and  temperatures.  The  influence  of  changing  bulk  strain  rate  on 
ignition  temperatures  was  much  smaller  than  the  role  of  turbulence  at  all  pressures. 
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The  observed  ignition  trend  can  be  explained  by  considering  discrete  mixing  by 
turbulent  eddies  similar  in  size  to  the  ignition  kernel,  which  can  bring  together  pockets  of 
reactants,  radicals,  and  high  temperature  gasses  that  would  not  otherwise  mix.  In  the 
second  limit,  where  ignition  is  governed  by  the  fast  reactions  that  define  the  crossover 
temperature,  slower  HO2  reactions  may  be  enhanced  by  the  re-introduction  of  these 
radicals  to  the  ignition  kernel.  This  new  source  of  chain  branching  would  provide  a  boost 
to  ignition.  As  the  turbulence  intensity,  and  thereby  frequency,  increases,  the  time 
available  for  this  new  branching  pathway  to  operate  decreases  and  the  ignition 
temperatures  increase  somewhat.  Similarly,  at  higher  pressures  where  ignition  is 
controlled  by  the  slower  HO2  pathways,  the  reduced  residence  time  due  to  turbulent 
fluctuations  will  hinder  ignition  such  that  turbulent  ignition  temperatures  become  higher 
than  in  laminar  flows. 

This  work  is  reported  in  Publication  No.  1 . 

1.3  Modeling  study  of  turbulent  nonpremixed  ignition 

The  above  experimental  study  on  the  effects  of  turbulence  on  nonpremixed  hydrogen 
ignition  showed  that  optimal  conditions  for  ignition  exist  in  weakly  turbulent  flows  where 
the  ignition  temperature  is  lower  than  either  laminar  or  strongly  turbulent  flows  at  the 
same  pressure.  In  order  to  imderstand  this  phenomenon,  computational  simulation  was 
conducted  using  Monte  Carlo  techniques  to  solve  a  joint  scalar  PDF  equation.  The 
simulation  involved  calculating  the  velocity  field  using  a  Reynolds  stress  model,  and 
using  this  frozen  flow  field  in  the  subsequent  ignition  calculation.  The  decoupling  is 
justified  because  the  small  amount  of  heat  release  prior  to  ignition  is  not  expected  to 
affect  the  flow  field.  Results  show  that  under  most  flow  conditions,  there  is  a  large 
increase  in  H  at  ignition  that  is  orders  of  magnitude  larger  than  the  background 
fluctuations.  Heat  release  is  shown  to  be  negligible  prior  to  ignition  by  examining  the 
stochastic  particle  temperatures  conditioned  on  the  mixture  fraction.  After  ignition  has 
occurred,  effects  of  heat  release  are  evident  in  the  conditional  temperature  data  and  the 
partially  decoupled  reactions  are  no  longer  expected  to  accurately  represent  the  flow 
conditions. 

Figure  3  shows  that,  at  atmospheric  pressure,  the  calculated  ignition  temperatures 
are  about  30K  higher  than  the  calculated  laminar  ignition  temperatures,  and  no 
appreciable  variation  is  noted  for  flows  with  different  turbulence  intensities.  Similar  over¬ 
predictions  were  noted  between  hirbulent  and  laminar  calculations  at  pressures  up  to  3-4 
atm.  At  still  higher  pressures,  the  turbulent  calculations  do  not  clearly  exhibit  the  third 
limit  behavior  found  in  laminar  ignition.  These  results  differ  from  the  experiment  which 
did  find  an  effect  of  turbulence  intensity  on  ignition  temperatures  and  show  a  decrease  in 
ignition  temperature  with  bulk  strain  rate  instead  of  the  increase  shown  in  the  turbulent 
calculations.  At  higher  pressures,  the  quantitative  agreement  between  turbulent 
calculations  and  experiments  improved  but  it  is  unclear  whether  the  trends  of  the  two  sets 
of  data  are  qualitatively  similar.  Further  modeling  work  is  being  conducted  aiming  to 
resolve  this  discrepancy. 

This  work  is  reported  in  Publication  No.  2. 
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1.4  Nonpremixed  ignition  ofH2/air  in  a  mixing  layer  with  a  vortex 
To  provide  further  insight  into  turbulent  ignition  in  general  and  optimize  ignition  in 
particular,  nonpremixed  ignition  in  an  elemental  structure  of  turbulent  flows,  namely  the 
vortex,  was  computationally  investigated  for  hydrogen/air  in  a  mixing  layer,  with  the  fuel 
side  being  pure  hydrogen  at  300K,  and  the  oxidize  side  being  hot  air  at  2000K. 

Figure  4  shows  the  maximum  heat  release  rate  for  mixing  layers  with  and  without 
(ID)  the  vortex.  It  is  seen  that,  compared  to  the  ID  case,  the  vortex  case  has  higher  heat 
release  rates  and  an  additional  peak,  hence  signifying  the  presence  of  two-staged  ignition. 
The  first  ignition  occurs  in  a  diffusion  mode  with  chain  branching  reactions  dominating. 
The  second  ignition  takes  place  in  the  premixed  mode,  with  more  chemical  reactions 
involved,  and  is  substantially  affected  by  the  heat  and  species  generated  in  the  first 
ignition  event.  Results  further  show  that  the  scalar  dissipation  rate  is  a  crucial  controlling 
factor  for  ignition,  which  causes  a  nonmonotonic  dependence  of  ignition  delay  on  all  the 
vortex  parameters. 

This  work  is  reported  in  Publication  No.  3. 


2.  Chemical  Kinetics  of  Laminar  Counterflow  Ignition 
2.1  Premixed  hydrogen/air  and  propane/air  ignition 

Nonpremixed  ignition  of  H2/Air  has  been  extensively  studied  in  previous  years,  yielding 
results  that  can  be  largely  interpreted  on  the  basis  of  homogeneous  systems.  While  these 
results  are  fundamentally  interesting  and  technically  important,  there  exist  many  practical 
situations  in  which  the  ignition  of  a  premixture  is  of  relevance.  Examples  include  the 
initiation  of  accidental  explosions  of  a  premixed  combustible  and  ignition  in  the  mixing 
layer  of  the  hot  re-circulation  zone  for  flame  holding.  Furthermore,  within  turbulent 
flames  there  are  pockets  of  rmbumed  combustibles  that  could  be  ignited  by  the  hot 
combustion  products.  Since  the  structure  and  response  of  nonpremixed  and  premixed 
flames  are  qualitatively  different,  it  is  therefore  of  interest  to  extend  our  studies  on 
nonpremixed  ignition  to  premixed  ignition. 

The  problem  studied  is  that  of  the  ignition  of  a  premixed  hydrogen/air  jet  by  a 
heated  counterflowing  inert  jet.  Figure  5  shows  that  the  ignition  temperatures  decrease 
when  the  mixture  equivalence  ratio  ((|))  is  reduced  from  0.5  to  0.2.  This  is  surprising  since 
the  mixture  reactivity  is  expected  to  be  higher,  and  the  ignition  temperature  lower,  as  (j) 
approaches  unity  for  the  present  lean  mixtures.  This  interesting  behavior  can  be  explained 
by  realizing  that  the  global  (|),  evaluated  at  the  cold  boundary,  does  not  directly  indicate 
the  mixture  stoichiometry  in  the  ignition  kernel,  and  hence  is  not  the  appropriate 
parameter  to  represent  the  system  stoichiometry.  The  ignition  kernel  is  located  on  the  hot 
side  of  the  flow  and  is  formed  by  the  diffusion  of  fuel  and  oxidizer  from  the  reactant  side. 
The  large  diffusivity  difference  between  Ha  and  O2  causes  the  more  diffusive  H2  being 
relatively  enriched  as  the  mixture  diffuses  towards  the  hot  flow.  As  a  result,  the  ignition 
kernel  is  much  richer  than  the  initial  reactant  mixture.  When  the  boundary  ^  becomes 
leaner,  the  local  <}»  in  the  kernel  moves  towards  stoichiometric,  the  mixture  reactivity 
increases,  and  the  ignition  temperature  decreases. 
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One  of  the  consequences  of  the  existence  of  the  oxygen-deprived  kernel  is  that 
reaction  rates  decrease  to  the  point  where  diffusive  transport  of  the  H  atom  from  the 
kernel  becomes  important.  The  diffusive  loss  is  most  pronounced  on  the  hot  side  of  the 
kernel  because  the  increasingly  rich  mixture  reduces  reaction  rates  further,  resulting  in  a 
tail  to  the  H  concentration  profile  towards  the  hot  boundary.  In  contrast  to  nonpremixed 
ignition,  the  diffusive  transport  is  similar  in  magnitude  to  the  chain  terminating  reactions, 
so  premixed  ignition  temperatures  are  typically  80K  higher  than  nonpremixed  ignition 
temperatures  and  are  more  sensitive  to  changes  in  the  strain  rate. 

The  dearth  of  oxidizer  in  the  ignition  kernel  can  be  remedied  by  adding  small 
amounts  of  oxygen  to  the  hot  nitrogen  flow.  It  was  found  that  the  ignition  temperatures 
approach  nonpremixed  values  by  adding  as  little  as  1%  O2  to  N2,  while  further  oxygen 
addition  has  little  effect  as  shown  in  Fig.  6.  Thus  by  increasing  the  reaction  rates  and 
reducing  the  diffusive  loss,  the  stoichiometry  of  the  premixed  kernel  with  oxygen 
addition  becomes  similar  to  that  of  the  nonpremixed  kernel  even  though  the  total  amoimt 
of  oxidizer  is  different.  This  therefore  constitutes  the  transition  between  premixed  and 
nonpremixed  ignition.  In  practical  situations,  it  is  reasonable  to  expect  that  some  amount 
of  the  oxidizing  species  would  be  present  in  the  hot  flow,  perhaps  as  a  result  of  lean 
combustion  or  the  dissociation  of  combustion  products,  and  the  scenario  of  premixed 
ignition  with  oxidizer  doping  is  a  possibility  that  needs  to  be  considered. 

This  study  was  fhrther  extended  in  two  directions:  hydrogen  ignition  at  elevated 
pressures,  and  propane  ignition  at  atmospheric  pressure. 

The  study  on  hydrogen/air  ignition  at  elevated  pressures  was  motivated  by  the 
recognition  that  ignition  phenomena  are  inherently  sensitive  to  detailed  chemistry,  and 
that  such  sensitivity  is  strongly  manifested  in  terms  of  pressure  variations.  Figure  7 
compares  the  characteristic  Z-shaped  pressure-temperature  explosion  limits  of  a 
homogeneous  mixture  with  the  present  premixed  counterflow  ignition  temperatures, 
determined  from  both  experiment  and  calculations,  for  9%  H2  in  air  mixture  ((|)  =  0.236) 
versus  heated  N2  at  a  pressure- weighted  strain  rate  of  k  =  216s'^  The  results  show  two 
important  distinctions  between  the  homogeneous  and  the  present  premixed  ignition 
situations.  First,  the  counterflow  ignition  response  curve  is  shifted  along  the  crossover 
temperature  path  to  a  higher  temperature  regime  due  to  the  reduced  residence  time. 
Second,  the  ignition  temperature  exhibits  five  limits  over  the  pressure  range  investigated. 
These  five  ignition  limits  have  the  following  characteristics. 

The  first  three  ignition  limits  qualitatively  mimic  the  homogenous  case. 
Quantitatively,  the  pronounced  H  radical  loss  to  the  hot  inert  side  pushes  the  first  limit  to 
end  at  a  higher  pressure.  The  steepening  of  the  second  limit,  which  turns  and  forms  the 
third  limit,  is  due  to  the  increasing  importance  of  R-18  (HO2  +  H2^  H2O2  +  H).  This 
reaction  converts  the  relatively  inert  HO2  to  the  more  reactive  H  radical  and  is  greatly 
facilitated  by  the  abundant  H2  within  the  ignition  kernel  owing  to  its  preferential 
diffusion  over  O2.  As  the  concentration  of  H  is  built  up  with  pressure,  reactions  between 
HO2  and  H  (RIO:  HO2  +  H2  +  O2  and  R1 1 :  HO2  +  OH  -i-  OH)  start  to  influence 
ignition.  The  combined  effect  is  to  weaken  the  termination  reaction  R9 
(H+02+M-»H02+M)  by  a  factor  of  2kl0/(kl0+kll).  Ignition  occurs  at  a  lower 
temperature  than  the  crossover  temperature,  and  an  extended  second  limit  behavior 
appears.  Finally,  as  more  HO2  are  produced  with  increasing  pressure,  its  recombination 
reaction  R14  (HO2  +  H02^  H2O2  +  O2)  generates  H2O2  which  further  decomposes  to 
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two  OH  through  R16  (H2O2  +  20H  +  M).  Thus  ignition  again  becomes  easier  at  the 

fifth  limit.  As  a  whole,  the  five  limits  of  premixed  ignition  can  be  interpreted  as  a  global 
three  limits  behavior,  with  the  third  (explosion)  limit  interrupted  by  the  extended  second 
limit  which  is  the  observed  fourth  ignition  limit. 

The  counterpart  of  hydrogen/air  ignition  is  propane/air  ignition  because  propane, 
being  less  mobile  than  oxygen,  will  render  an  originally  lean  mixture  more  fuel  lean  in 
the  ignition  kernel.  Thus  while  increasing  (j)  will  facilitate  ignition,  the  extent  is  expected 
to  be  weakened.  Figure  8  shows  that  the  ignition  temperature  decreases  with  increasing  (|), 
although  the  influence  of  preferential  diffusion  cannot  be  readily  separated  out. 

This  work  is  reported  in  Publication  Nos.  4  and  5. 

2.2  Counterflow  ignition  temperatures  and  laminar  flame  speeds  of  C2-C3  hydrocarbons 
at  atmospheric  and  elevated  pressures 

While  extensive  kinetics  data  exist  for  hydrogen  and  methane,  comparatively  fewer  data 
have  been  acquired  for  the  hierarchically  more  complex  C2  and  C3  hydrocarbons,  which 
are  important  fuels  and  intermediate  constituents  in  hydrocarbon  combustion.  Therefore, 
the  primary  objective  of  the  present  study  is  to  acquire  additional  data  on  the 
nonpremixed  counterflow  ignition  temperatures  as  well  as  the  laminar  flame  speeds  of 
most  of  the  C2-C3  hydrocarbons,  namely  acetylene  (C2H2),  ethylene  (C2H4),  ethane 
(C2H6),  propylene  (C3H6)  and  propane  (C3H8).  These  data  cover  the  regimes  of  high- 
temperature  oxidation  in  flames  and  low-  to  intermediate-temperature  oxidation  relevant 
for  ignition,  in  premixed  as  well  as  nonpremixed  situations,  and  with  substantial  pressure 
variations.  These  data  were  then  compared  with  those  calculated  by  using  an  optimized 
mechanism  for  C1-C3  (Z.W.  Qin,  V.V.  Lissianski,  H.X.  Yang,  W.C.  Gardiner,  S.G.  Davis 
and  H.  Wang,  Proc.  Combust.  Inst.  28,  1663-1669,  2000)  and  an  ethylene  ignition 
mechanism  (B.  Varatharajan  and  F.A.  Williams,  J  Propul.  Power  18,  344-351,  2002). 
These  two  mechanisms  are  referred  to  as  the  Q-  and  VW-mechanisms  respectively. 

Results  (Figs.  9  and  10)  show  that,  for  all  ignition  cases,  the  Q-mechanism 
consistently  predicts  higher  ignition  temperatures  than  the  experiments,  while  the  VW- 
mechanism  agrees  well  with  the  experiments.  On  the  other  hand,  the  Q-mechanism 
agrees  with  the  laminar  flame  speed  data  moderately  well.  For  acetylene  and  ethylene, 
the  Q-mechanism  slightly  over-predicts  and  under-predicts  the  lean  and  rich  mixtures, 
respectively,  as  shown  in  Fig.  1 1  and  12.  A  possible  source  for  the  stronger  reactivity  of 
the  Q-mechanism  is  that  it  was  optimized  against  previous  ethylene  flame  speed  data 
based  on  linearly  extrapolation,  which  yields  higher  laminar  flame  speeds.  A  fairly 
detailed  exploration  was  conducted  on  the  causes  for  the  various  disagreements  between 
experiment  and  calculations. 

This  work  is  reported  in  Publication  No.  6. 

2. 3  Ignition  of  dimethyl  ether  in  counterflow 

Dimethyl  Ether  (DME),  CH3OCH3,  is  an  attractive  alternate  fuel  and  fuel  additive  for 
diesel  engines.  The  combination  of  its  high  cetane  number  and  low  boiling  point  of -25°C 
facilitates  mixing,  ignition,  and  cold  start.  The  pyrolysis  and  oxidation  process  of  DME 
have  been  studied  mainly  based  on  homogeneous  systems,  which  has  led  to  the 
development  of  two  detailed  mechanisms  reported  respectively  in  1998  and  2003. 
However,  their  validity  has  not  been  experimentally  scrutinized  in  nonhomogeneous 
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systems,  particularly  at  low  temperatiures  relevant  for  ignition.  In  the  present  study  the 
ignition  temperatures  of  DME  in  nonpremixed  counterflow  were  experimentally 
measured  and  subsequently  used  to  evaluate  the  accuracy  of  these  two  mechanisms. 

Specifically,  the  ignition  temperatures  of  nitrogen-diluted  DME  by  heated  air  in 
counterflow  were  experimentally  determined  for  DME  concentration  from  5.9  to  30%, 
system  pressure  from  1.5  to  3.0  atmospheres,  and  pressure-weighted  strain  rate  from  110 
to  170s'^  Figure  13  shows  representative  results  for  the  ignition  temperature  as  a 
function  of  DME  mole  firaction  in  nitrogen  at  3  atm  and  a  pressure- weighted  strain  rate  of 
150s'\  It  is  seen  that  the  ignition  temperature  monotonically  decreases  with  increasing 
fuel  concentration  in  the  fuel  jet,  and  that  such  a  dependence  becomes  stronger  for  lower 
fuel  concentrations.  Furthermore,  while  the  1998-mechanism  substantially  over-predicts 
the  ignition  temperature,  the  2003-mechanism  agrees  very  closely  with  the  experimental 
data.  This  observation  holds  for  all  results  obtained  in  the  present  study. 

Sensitivity  analysis  was  performed  to  identify  the  differences  between  these  two 
mechanisms.  Overall,  the  1998-mechanism  shows  much  reduced  sensitivity  to  the  low 
temperature  oxidation  sub-mechanism,  while  the  2003-mechanism  shows  that  the 
majority  of  the  sensitivity  comes  firom  the  low-temperature  kinetics.  For  example,  the 
2003-mechanism  exhibits  a  much  stronger  reactivity  around  600K,  as  manifested  by  the 
turning  point  behavior  (Fig.  14).  On  the  contrary,  the  1998-mechanism  indicates  the  CH3 
radical  concentration  grows  with  the  hot  boundary  temperature  only  gradually,  which 
produces  considerably  less  low  temperature  reactivity  and  therefore  causes  the  ignition 
temperature  to  be  1039  K,  which  is  about  100  K  higher  than  that  predicted  by  the  2003- 
mechanism.  The  reason  that  the  low  temperature  kinetics  contributes  so  much  to  ignition 
is  that  before  the  fuel  enters  the  high-temperature  ignition  region  around  the  hot-stream 
boundary,  it  is  progressively  heated  from  its  low-temperature  region. 

The  2003-mechanism,  consisting  of  79  species  and  398  elementary  reactions,  was 
then  systematically  simplified  by  using  directed  relation  graph  theory  to  a  skeletal 
mechanism  of  49  species  and  251  elementary  reactions,  which  in  turn  was  further 
simplified  by  using  computational  singular  perturbation  and  quasi-steady-state  species 
assiunption  to  a  reduced  mechanism  consisting  of  33  species  and  28  lumped  reactions.  It 
was  demonstrated  that  both  the  skeletal  and  reduced  mechanisms  mimic  the  performance 
of  the  detailed  mechanism  with  high  accuracy,  as  shown  in  Fig.  15. 

This  work  is  reported  in  Publication  No.  7. 
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Figure  1 :  Effect  of  turbulence  intensity 
on  ignition  temperature  at  latm 
and  20%  H2  concentration. 
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Figure  3;  Comparison  of  the  effects  of 
different  mixing  models  on  the 
ignition  response  to  changing 
pressure. 


Figure  2:  Effect  of  pressure  on  ignition 
temperature  for  different 
tvurbulence  intensities. 
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Figure  4:  Effect  of  turbulence  intensity 
on  ignition  temperature  at  latm 
and  20%  Ha  concentration. 
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Figure  5:  Critical  temperatures  and  local 
(|)  as  functions  of  boundary 
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Figure  7:  Comparison  of  the  effects  of 
different  mixing  models  on  the 
ignition  response  to  changing 
pressure. 


Figure  6:  Effect  of  adding  O2  to  the  hot 
side  on  the  ignition 
temperature. 
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Figure  8:  Effect  of  turbulence  intensity 
on  ignition  temperature  at  latm 
and  20%  H2  concentration. 


15 


Local  Equivalence  Ratio 


Pressure  (atm) 


Figure  9:  Effect  of  pressure  on  ignition 
temperatures  for  C2H4  and 
C3H6.  Conditions:  5.9%  Fuel  in 
N2andk=150s‘‘ 
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Figure  1 1 ;  Laminar  flame  speed  of 

acetylene  (C2H2)  in  air  at  latm 
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Figure  10:  Effect  of  pressure  on  ignition 
temperatures  for  C2H6  and 
C3H8.  Conditions:  5.9%  Fuel 
inN2andk=150s-l 
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Figure  12:  Laminar  flame  speed  of 

ethylene  (C2H4)  in  air  at  latm. 
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Figure  13:  Effects  of  DME  mole  fraction 
on  ignition  temperature. 
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Figure  14:  Calculated  S-curves  for  30% 
DME  in  N2  mixture  with  150  s' 
'  strain  rate  at  3  atm,  using  both 
mechanisms 


Figure  15:  Comparison  of  the  computed 
results  using  detailed,  skeletal 
and  reduced  mechanisms  for 
the  ignition  temperatures. 


17 


